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Abstract 

We study the radiative electroweak symmetry breaking and the relic abundance of neutralino 
dark matter in supersymmetric type I seesaw model. In this model, there exist threshold cor- 
rections to Higgs bilinear terms coming from heavy singlet sneutrino loops, which make the soft 
supersymmetry breaking (SSB) mass for up-type Higgs shift at the seesaw scale and thus a min- 
imization condition for the Higgs potential is affected. We show that the required fine-tuning 
between the Higgsino mass parameter ^ and SSB mass for up-type Higgs may be reduced at elec- 
troweak scale, due to the threshold corrections. We also present how the parameter /i depends 
on SSB B-parameter for heavy singlet sneutrinos. Since the property of neutralino dark matter is 
quite sensitive to the size of /x, we discuss how the relic abundance of neutralino dark matter is 
affected by the SSB B-parameter. Taking the SSB B-parameter of order of a few hundreds TeV, 
the required relic abundance of neutralino dark matter can be correctly achieved. In this case, dark 
matter is a mixture of bino and Higgsino, under the condition that gaugino masses are universal 
at the grand unification scale. 

PACS numbers: ll.30.Qc, ll.30.Pb, 12.60.Jv 
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I. INTRODUCTION 



Supersymmetric (SUSY) seesaw model is a SUSY extension of seesaw model [l, 2| which 
naturally explains small masses of neutrinos and stabilizes the hierarchy between electroweak 
scale and some other high scale without severe fine-tuning if the mass spectrum of super- 
partners are less than TeV scale as well. In SUSY type I seesaw model, we introduce not only 
heavy right-handed(RH) Majorana neutrinos but also their super partner called sneutrinos 
which are standard model gauge singlet. This leads us to anticipate that some predictions 
of MSSM can be deviated due to the contributions associated with the heavy RH neutrinos 
and their super partners, and new phenomena absent in MSSM may occur in SUSY type I 
seesaw model. With this regards, there have been attempts to study lepton flavor violation 
and neutrino masses in SUSY type I seesaw model 3-^. On the other hand, the gauge singlet 
RH neutrino superfield may affect Higgs sector as investigated in Ref. gI, where they have 
shown that there is a sizable negative loop contribution to the mass of the lightest Higgs 
field in the split-SUSY scenario at the price of giving up the naturalness in supersymmetry. 

In this study, we revisit the issue as to how the Higgs sector can be affected by heavy 
singlet sneutrinos while keeping naturalness in supersymmetry. It is well known that the 
lightest CP-even Higgs mass in the MSSM can get large one-loop corrections which increase 
with the top quark and squark masses THlQl]. The current experimental bound on the 



lightest C 
500 GeV 



-even Higgs mass, ruh ^114 GeV, demands top squark mass to be larger than 



111], which in turn leads to a fairly large correction to the soft supersymmetry 



breaking (SSB) mass for the up-type Higgs rajj^. In the MSSM, electroweak symmetry can 
be broken due to the large logarithmic correction to Tnjj^ 12|-|l6|. However, as is known, 
we need rather large fine-tuning between the Higgsino mass parameter fi and SSB mass 



m 



"jj^ to achieve the Z-boson mass at the electroweak scale through a minimization condition 
for the Higgs potential of the MSSM . In this study, we show that there exist some new 
contributions generated from the loops mediated by heavy singlet sneutrino sector to SSB 
mass mjj^ and the Higgsino mass parameter /i in SUSY type I seesaw model. The new 
contributions are given in terms of SSB parameters -Bat and SSB mass term for the singlet 
sneutrino at the seesaw scale. 

Integrating out the singlet neutrino superfield below the seesaw scale, SUSY type I seesaw 
becomes equivalent to the MSSM but those new contributions are taken to be threshold 
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corrections to Higgs bilinear terms. As will be discussed, those threshold corrections can 
lower the sizes of m^^^ and fi at the electroweak scale and thus the fine-tuning may be 
reduced. This means that the fine-tuning required for the radiative electroweak symmetry 
breaking can be shifted to tuning the size of -Bjv at the seesaw scale. In this paper, we 
investigate how the sizes of and fi at the electroweak scale depend on the parameter 

Since the property of neutralino dark matter is quite sensitive to the size of fi, we discuss 
how the relic abundance of neutralino dark matter is affected by the parameter B^. In 
fact, there exist some literatures in which the impacts of neutrino Yukawa couplings on 



neutralino dark matter in SUSY type I seesaw model have been discussed 17H24l|. It has 
been found that some regions of parameter space can significantly affect the neutralino 
relic density without the threshold corrections associated with the heavy singlet neutrino 
superfield. In our work, however, we consider possible existence of the threshold corrections 
generated from the loops mediated via the heavy singlet neutrino superfield which can also 
significantly affect the neutralino relic abundance by lowering the sizes of m^^^ and fi at the 
electroweak scale. Such a possibility of the impact on the neutralino relic density has not 
been studied before. 

This paper is organized as follows. First, we present the effective potential for Higgs 
fields in SUSY type I seesaw model in section II. We show that threshold corrections to 
Higgs bilinear terms are generated from the loops mediated by heavy singlet neutrino su- 
perfields. In section HI, we give the alternative derivation for the threshold corrections, 
using renormalization group equations (RGBs) for a general field theory. In section IV, we 
study the contributions of the threshold corrections to the radiative electroweak symmetry 
breaking and investigate how the size of the parameter fi can be affected by them. In section 
V, we discuss the relic abundance of neutralino dark matter. Finally section VI is devoted 
to conclusions and discussions. The details of convention for CP phases and derivation of 
the effective potential for Higgs fields are given in Appendix. 

II. THE EFFECTIVE POTENTIAL OF SUSY TYPE I SEESAW MODEL 

In this section, we first derive the effective potential of SUSY type I seesaw model, and 
then show that there exist threshold corrections to Higgs bilinear terms arisen due to the 
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heavy RH singlet sneutrinos. Those threshold corrections may be modified by wave function 
renormalization for Higgs field. 

The super potential of the SUSY seesaw model is given by 

W = fiH^-H2- Y,{L ■ H2W - ^iV^iV^ (1) 

where N'^ is a gauge singlet chiral superfield, which contains a RH neutrino and its scalar 
partner. Mr denotes the mass of the RH neutrino. Here, we do not consider the terms asso- 
ciated with the charged leptons and quarks whose contributions to our study are negligibly 
small except for top quark superfield. From now on, we consider only one generation of N'^ 
for simplicity, and the extension to three generations is straight-forward. The soft breaking 
terms of the Lagrangian in SUSY seesaw model are given by 

Aoft = ~ml\L\''-ml\N\''-{^B*^MlN'' + h.c. 

+2Re(5/ii7i ■ H2) - ml^HlHi - m]j^HlH2 

+ {a,Y,{H2 ■ L)N* + h.c^ , (2) 

where we can take Mr^Bn^Y^^h to be real by superfield rotation and U{1)r symmetry, 
whereas Ay and B are left as complex numbers. We discuss the details of the phase con- 
vention in Appendix A. From the superpotential given in Eq. ([T]), the SUSY part of the 
Lagrangian is obtained as follows: 

-Csusy = -KL ■ H2 + MnN*\^ - \YyN*L - fiH^l^ - \fi\^HlH2 - Y^\N\^hIH2 

- ]^MRN^Nn' - YJN^Il ■ + h.c. (3) 

With this Lagrangian, we can derive the effective potential by using field dependent 
masses for the singlet RH neutrinos and sneutrinos. The effective Higgs potential which 
includes 1-loop contributions mediated by the singlet RH neutrino superfields is written as 



+ (^log -^K + ml + \Al\) + 2ml + 2Re.(A,i?^) ) hIh^ 

- 2Re. \j^{Bn + a, log ■ H^j - £d, 



(4) 



where Q is a renormalization scale and Cd is D-term contributions given by 

^D = -Y (^i^^i - ^2^2)' - J [hIt^H, + Hlr^H.y . (5) 

In Appendix B, we present in detail how the effective potential is derived. Matching this 
effective potential with that of MSSM at the seesaw scale, we can obtain some relations 
between MSSM parameters and corresponding ones in SUSY seesaw model. Here, we do 
not include the loop contributions mediated by top quark and its super partner because 
they are identical to each other in both MSSM and SUSY seesaw model, and thus canceled 
in the relations. Therefore those contributions are irrelevant to the threshold corrections for 
the Higgs bilinear terms. The Higgs potential of the MSSM is given by, 

- (^B{Q^)fiHf ■ + h.c^ - Cd. (6) 

By matching the Higgs potentials Eq.(l6]) with Eq.(j4]) at = M^, we obtain the following 
relations, 

B{Ml) = B{Ml) + -1^5^. (7) 

On the other hand, the wave function renormalization for the Higgs field H2 in the limit of 
small external momenta is given by 

{''-^^'°^'§)''''"^^''"?' 

where we neglect the terms suppressed by M^^. We notice that there exist no contributions 
from heavy RH neutrino superfields to wave function renormalization for Hi. At = Af^, 
Eq.® becomes dfj,H\d^H2, so the relations given in Eq.d?]) are not modified by wave function 
renormalization. 

It is worth noting that the soft breaking parameter of singlet sneutrino, -Bjv, contributes 
to the Higgs mass fn\j^{M'j^ and the parameter B. We use RGEs for the soft breaking 
parameters of the MSSM to obtain their low energy values below the seesaw scale Mr, 
whereas the corresponding RGEs given in the SUSY seesaw model are used above the seesaw 
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scale. Thus, the values of the parameters in the RH side of Eq.([7]), mjj_^{Q'^ = M^) and 
m\^{Q'^ = M|,), depend on the boundary condition at further high energy scale, such as 

MgUT or Mpianck- 



III. THE THRESHOLD CORRECTIONS FROM RENORMALIZATION GROUP 
EQUATIONS 

In this section, we study the alternative derivation of the threshold corrections given in 
Eq.(I7j) by using RGEs including threshold effects. The RGEs in MSSM including threshold 



effects are discussed in Refs. 
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-|28|. We derive the one-loop RGEs for Higgs mass squared 



parameters in the SUSY seesaw model, by using the formulas for RGEs of dimensional 



parameters in general gauge field theories 29|. Then we integrate them and obtain the 



threshold corrections. Here we focus on the effects from the heavy neutrino and sneutrinos. 

The key point of the derivation of the threshold corrections is to take into account three 
different thresholds. One of them corresponds to the mass of RH neutrino (Mr), and the 
others correspond to the masses of the heavy sneutrinos , i.e., the super partners of the RH 
neutrino. They are two real scalar fields and their masses are deviated from due to soft 
SUSY breaking terms of the sneutrinos sector , as given by 

Cra^. = -ImIn', - ^MlNl (9) 

where Ni and N2 are real and imaginary part of the complex scalar field A^, respectively 
and are defined as, 

Ni = {N + N*)/V2, N2 = {N -N*)/{V2i). (10) 
The masses of the A^i and N2 are then give by 

= + + BnMr, Ml^ = m% + Ml- BnMh. (11) 

Since B]\f is real positive, the hierarchy of the three mass scales is given by 

Ml >Ml> Ml^. (12) 

Then the energy scales at which A^^i, N2 and are decoupled are different each other, 
yielding the threshold corrections to Higgs mass squared parameters. The Higgs mass terms 
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are given as 



Higgs 



(13) 



where, 



mil = 


1 |2 1 2 


"^22 = 


1 |2 1 2 

|/i| 


m?2 = 





(14) 



Following 29|, we divide all the complex scalar fields into their real and imaginary parts, 
and derive the beta functions for the Higgs mass squared parameters by adopting the step 
functions of the renormalization scale {Q) to take into account the thresholds. Then we 
obtain the threshold corrections by integrating the beta functions with respect to the energy 
scale between two mass scales of the singlet sneutrinos. 

At one-loop level, the beta functions for the Higgs mass parameters are given as. 



(4vr: 
(An 



(4vr: 



.2 dml^ 
d\nQ 

2 dmT^2 



2 ^^22 

ding 



Ml 



N2' 



Y^'A^fj. 9{Q'-Ml) + 9{Q'-Ml: 



-Y^'fiMn9{Ml^-Qy{Q 



N2^ 



+ e{Q^ - M' 



N2' 



-Y^ [2Re(A,) + Bn] Mn OiM^ - Q^iQ' - M^] 



e{Q^-Ml) + e{Q'-Ml^)-2e{Q' 

+2r>2^ 9{Q' - Ml) + Ylm\ \e{Q^ - M|J + e{Q^ 



Ml) 



Mi 



N2' 



(15) 



Here, we note that only the terms coming from the neutrino-sneutrino sector are presented 
because the other terms are the same as those in MSSM. In deriving the RGEs, we take into 
account the fact that the effective theory changes as passing each threshold corresponding 
to the heavy degree of freedom. At the energy scale above where the RH neutrino and 
sneutrinos are active, our RGEs given in Eg. (1151) are consistent with those in supersymmetric 



type I seesaw model 
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3l| . While the RH neutrino and the lighter sneutrino are active 



between the two scales Mjy^ and Mr-, only the lighter sneutrino is active between the two 



scales Mr and M^^. Finally, the effective theory becomes MSSM below M^^. In each step. 



N2- 



we integrate out the heavier degrees of the freedom and derive the effective theories which 
are vahd at the lower energy scales. 

By integrating the beta functions with respect to Q from down to M^^, we obtain 
the threshold corrections. Since the integrals can be approximated as follows; 

Bn 



cilnQ 



In 



M 



N2 

M 



Mm, 



JVi 



d\nQ 



In^^ 

Mn 



Mf 

1 
2 



B 



N _|_ nT'N2 



Ml 



2M| 



(16) 



Only the terms proportional to Mji or in Eq.( IT5l) contribute to the threshold corrections. 
The results of integrating the beta functions give 



6B 



Y 



^^[m% + ReiA,)B^]+0{M],'] 



(17) 



which are the same as Eq.Q. 

Next, we discuss how the numerical value of the parameter yU can be affected by threshold 



corrections 
scenario 
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br the Higgs bilinear terms in the radiative electroweak symmetry breaking 



16| . In the calculation, we assume that gaugino masses, scalar masses and 



A-terms are universal at the GUT scale. 



IV. MU TERM AND RADIATIVE ELECTROWEAK SYMMETRY BREAKING 

As we have shown, the soft breaking parameter for the Higgs mass m^^^ in the MSSM at 
the seesaw scale Mr is determined by not only fn\j^{M'j^ calculated via RGEs in the SUSY 
seesaw model but also additional contribution due to the loops mediated by light and heavy 
sneutrinos in the seesaw model at the scale M/j . From Eq.Q, the shift of Tn\j^ from fh\^ 
at the scale Mr is approximately given as. 

Therefore the soft breaking parameter Bn of the order of 500 TeV may significantly affect 
at the scale Mr. This observation in turn indicates that the shift of m'jj^ at the scale 
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Mr affects electrowealc symmetry brealfing in the MSSM when we iske the MSSM as an 
effective theory of SUSY type I seesaw model at low energy scale. 

Let us discuss how electroweak symmetry breaking can be affected by the parameter B^q. 
In the MSSM, radiative breaking of electroweak symmetry can occur when SSB parameters 
for Higgs sectors satisfy the following relation: 



2-z = -1/^1+ • (19) 



In the limit of large tan /3, this relation becomes 

1 o 



2^. 



/i|^-m^,(mi). (20) 



Therefore we see that the value of /i and directly related. In order to satisfy this 

condition, has to be negative at the scale irtz- In the radiative electroweak symmetry 
breaking scenario, m|^^ is generally taken to be positive at high energy scale, but it receives 
quite large radiative corrections due to heavy stop mass and large top quark Yukawa cou- 
plings between high and low energy scales, which drive m^^^ negative so that electroweak 
symmetry can break at low energy scale. At the scale above Mr, soft breaking masses and 
couplings are subject to the RGEs of the SUSY seesaw model. The RGE for Tn\j^ in the 



SUSY seesaw model is given by [30|, ISl] 



dmij^ _ 2 



dt 167r 



2 



-\glMl - ?>glMl + 3r/Xi + Y^X„ 
5 



(21) 



where t = In = rrA + m? + rni,^ + lA^P and X„ = + + rrii,^ + \AA'^. Here, 

VO V3 *-R 112 ' ' L N 112 I I 

Ml and M2 denote the bino mass and the wino mass, respectively. The last term comes 
from the presence of RH neutrino superfields and other terms are the same as those in the 
MSSM. 

It is expected that the RGE for mjj^ can be significantly affected by the Yukawa coupling 
of neutrino sector when it is quite large. We can estimate the deviation of the from 
that without neutrino sector by integrating out eq. (12T]) explicitly. The deviation at the scale 
Mr is approximately given as, 

Sio,ml^^^^{3ml + Al)lnj^, (22) 

where mo and Aq are the universal values for scalar masses and A-terms respectively. For 
Mr = 6 X 10^'^ GeV and Mx ~ 2 x 10^^ GeV, this contribution can be written approximately 



5w< - -5.5 X lO^(GeV)^ (^)^ (23) 

As we can see from eq. (fT8l) . Siogfn'j^^ is easily dominated by the threshold correction when 
-Bat is large. 

Without threshold corrections, the weak scale value of m^^^ becomes more negative than 
that of minimal supergravity (mSUGRA) case. This affects the condition for electroweak 



symmetry breaking and the allowed regions for the observed relic density of dark matter 



13, 



18| . Especially, the allowed region where is small, is changed significantly. Universal 



scalar mass at the GUT scale, mo is larger than that of mSUGRA. However with inclusion 
of the threshold corrections, mo can be smaller than that of mSUGRA when is large. 

Figure [H shows the RG evolution of m^^ m? with the energy scale. Here, m^- is 
defined as m? = m^m^-^. We assume that soft breaking masses, gaugino masses and A- 
terms are universal at the GUT scale(~ 2 x lO^^GeV). The calculations are performed with 
ISASUGRA code which is included in ISA JET package 32]. The input values used in the 
calculations are given in the caption and neutrino masses rrii, and M/j are taken to be 0.1 eV 
and 6 x 10^'^GeV, respectively, in both panels so that Y,^ and Yf become the same order of 
magnitude. The pink, blue and red curves correspond to the predictions of sign{ml[^)\mH^\ 
including threshold corrections for Bn = 500, 50 and 5 TeV, respectively. The green curves 
show how the predictions of m? evolve from the GUT scale to the electroweak scale. When 
Bn = 50 TeV, A,y ~ 300 GeV and mo ~ 1 TeV, the threshold correction and the running 
effects from the neutrino Yukawa sector are almost canceled, i.e. Sm'jj^ + ^log^'u^ ~ 0- 
Therefore the blue lines below the scale behave as if there are no effects from neutrino 
Yukawa sector. As we can see from Fig. [H the value of m|^^ at the scale m^ obtained in 
the SUSY seesaw model is significantly deviated from that obtained in the MSSM for given 
input values of mo, mi/2, ^o, tan/3 and Bj^ = 500TeV, whereas such a deviation disappears 
for Bn < 5TeV. 

In the case without threshold corrections, the running of the mj^ „ in mSUGRA with RH 
neutrino superfield (mSUGRA+RHN) is discussed in Refs. 



18 



23H. The weak scale values 



of yln^nj tend to be larger than those in mSUGRA scenario. The difference between 
mSUGRA and mSUGRA+RHN is up to a few hundred GeV, when mo > 1.5 TeV and 
Yiy > Yt. On the other hand, our results show that the threshold correction increases 
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mH^{Q'^ = M^) by several hundred GeV and therefore the weak scale values of can 
be smaller than those in mSUGRA scenario when Bf^r is large. 

The significant deviation of mj^^ at the scale mz in turn leads to a significant change in 
|yu| through the stationary condition, Eq. flTI?]) . In Fig. [21 we present how \fi{Mz)\ depends 
on the value oi B^. As the value of B^ increases, becomes smaller, due to the threshold 
corrections to rn?fj^{MR). 

It is worthwhile to notice that the size of the mass parameter /i characterizes the property 
of neutralino dark matter. Since /i is the Higgsino mass term, changing /i may affect the 
composition of the neutralino dark matter. This indicates that relic abundance of the dark 
matter is affected by -Bat, especially on the condition that gaugino masses are universal at 
the GUT scale. 



V. BINO-HIGGSINO DARK MATTER 



In this section, we show that the lightest SUSY particle (LSP) is a bino-Higgsino mixture 
state when the size of parameter Bn is of the order of several hundred TeV, and the result of 
the WMAP observation can be accounted for well. Here, we assume that soft scalar masses, 
gaugino masses and A terms are universal at the GUT scale. We consider the lightest 
neutralino dark matter candidate. 

The neutralinos are the physical states which are composed of the bino, wino and two 
Higgsinos. The neutralino mass matrix in the B—W—H1—H2 basis is given by, 

/ Ml —mz cos (3 sin 6w mz sin {3 sin 6w \ 

M2 mz cos /3 cos 6y/ —mz sin (3 cos 6w 

-mz cos (3 sin 6y/ mz cos (3 cos 6w — /i 

y sin /3 sin 6*1^ — sin /3 cos 6'vi/ —[^ J 

where Mi and M2 are the bino and wino masses, respectively, and 9w is the Weinberg angle. 
This matrix is diagonalized by the unitary matrix A^, 

j^d^ag ^ N*M,,N-\ (25) 

In terms of A^, the lightest neutralino is expressed as a mixture of the gauginos and the 
Higgsinos: 

,0 



,(24) 



X 



NnB + N12W + A^iaifi + NuH2. 



(26) 
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FIG. 1: The renormalization group evolutions of soft scalar masses for up-type Higgs 
and stops are shown. The calculation is performed by taking mp, ^O) tan /? to be 
ITeV, 400GeV, SOOGeV, 10, respectively in the upper panel and 700GeV, 500GeV, SOOGeV, 20, re- 
spectively in the lower panel. We take neutrino masses and M/j to be O.leV and 6 x 10^'^GeV, 
respectively, in both figures so that Yy and Yt become the same order of magnitude. The pink, 
blue and red curves correspond to the predictions of sign(m\i^\mH^\ for Bj^ = 500,50 and 5 TeV, 
respectively. The green curves correspond to the MSSM prediction of m^. 
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FIG. 2: The values of are plotted as a function of The lower red line is obtained for 

mo = ITeY ,1711/2 — 400GeV,^o = 300GeV and tan/3 = 10, and the upper green line for ttt-q = 
700GeV,mi/2 = 500GeV,^o = 300GeV and tan/3 = 20. We take the same values of m^, and Mr 
as in Fig. 1. 

Since we assumed a universal value for gaugino masses at the GUT scale, gaugino masses 
Mi are related to gauge couplings Qi as follows; 

M.(Q) gKQ) 



M{Agut) 9^{Agut)' 
and this relation is easily derived from renormalization group equations for gauginos. 



(27) 



(28) 



dt 167r2"^^' 

where 6j are coefficients of beta-functions for Qi. From Eq. fl27|) . the bino mass Mi is written 
in terms of the wino mass M2: 

5 



Ml = - tan^ ewM2 ^ O.SMs, 
3 



(29) 



at the scale m^. 

The relic density of a cold dark matter, f2cDM^^, is determined by WMAP observation 



33| and its value is given by 



QcDivi/i^ = 0.1131 ±0.0034. 



(30) 
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For |/i| 3> M2, the dark matter is bino-like, whereas for \fi\ ^ M2 the dark matter is 
Higgsino-hke . In general, a bino-hke dark matter leads to a large relic abundance of a dark 
matter, which can not accommodate the result from WMAP observation. This is because 
couplings for bino are smaller than those for Higgsino and wino. When the value of 
decreases, the Higgsino fraction defined by jA'^isP + |A'^i4p increases, which leads to larger 
annihilation cross sections for Higgsino-like dark matter. Therefore we can fit the right 
amount of relic abundance derived from the result of WMAP observation with dark matter 
candidate composed of a bino-Higgsino mixture. 

As we can see from Eq.(|T9l). since the value of |/i(m^)p becomes smaller as Bn{Mr) 
increases. Larger value of Bj:q{M^) leads to larger Higgsino fraction, which makes the relic 
abundance of a dark matter decreased. Fig. [3] presents the predictions of relic abundance 
of the lightest neutralino and corresponding contributions of Higgsino components as a 
unction of B^- Our numerical calculation is performed by using micrOMEGAs 2.2 code 
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351]. The blue line represents the value of the relic abundance obtained from WMAP 



observation. In this figure, we can see that as Bn{Mii) increases, Higgsino fractions get 
larger, which makes relic abundances smaller. From our numerical analysis, it turned out 
that the right amount of the relic abundance of the dark matter could be explained by 
taking the parameter B^ to be of the order of several hundred TeV which makes Higgsino 
fractions large. The allowed regions of parameter space for the observed relic density of 
the dark matter are most conveniently shown in (mi/2,mo) plane. In mSUGRA+RHN 
scenario without threshold correction, the allowed regions are given in Refs. 18|, l23| . One 
of the regions corresponding to the small is located along the region where electroweak 
symmetry breaking can not take place. This region corresponds to mg > 1.3TeV. The 
values of mo depend on the renormalization group running effect from neutrino Yukawa 
sector and it decreases the low energy value of Tn\j^. When this effect becomes larger, we 
need to choose larger mo as the GUT boundary condition. With inclusion of the threshold 
correction to mf^^; however, the consequences change. In our scenario, as shown in Figj3l 
we can take mo as small as 700GeV, since the threshold correction is added to m^^^ at the 
scale Mr. Therefore, we conclude that the allowed regions where the observed relic density 
is explained by the bino-Higgsino dark matter are very different from those of mSUGRA 
and mSUGRA+RHN scenario. 
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FIG. 3: The relic abundances of the Hghtest neutrahno (red curves) and correspond- 
ing Higgsino contributions (green curves) are drawn as a function of B^. We take 
"^0l 'TT'i/2i ^0) tan/3 to be ITeV, 400GeV, SOOGeV, 10, respectively in the upper panel and to be 
700GeV, 500GeV, 300GeV, 20, respectively in the lower panel, fi is positive, and the values of nii, 
and Mfi are taken to be the same as in Fig. 1. The blue lines correspond to the value of the relic 
abundance obtained from WMAP observation. 
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VI. CONCLUSION AND DISCUSSION 



We have investigated the effective low energy Higgs potential of the SUSY type I seesaw- 
model. We found that Higgs bihnear terms got threshold corrections at the scale below M^j, 
due to heavy singlet sneutrino loops. These threshold corrections are proportional to B- 
term of heavy singlet sneutrino Bn- Therefore if is large enough, the mass parameters 
of Higgs bilinear terms arc significantly shifted at the scale M^, which in turn leads to 
a shift of the parameter and reduction of the fine-tuning between the Higgsino mass 
parameter n and SSB mass for up-type Higgs at the electroweak scale. We presented how 
the parameter /i depends on Bj^. We have shown that dark matter becomes a mixture 
of bino and Higgsino for B^ oi the order of several hundreds TeV and the observed rehc 
abundance can be consistently explained by the bino-Higgsino dark matter. It turned out 
that the allowed region of parameter space constrained by the relic abundance of dark matter 
in this model is very different from the MSSM without seesaw under the assumption that 
SSB terms are universal at the GUT scale, mainly because of the threshold corrections to 
mjj^ . Our results are also different from those of conventional mSUGRA with type I seesaw 
which does not include the threshold corrections to rnir . 

Naturalness problem for such a large value of Bn is beyond the scope of this work. 
Since the size of Bn of order of several hundreds TeV is much larger than the scale of soft 
breaking parameters, the origin of B^ must be different from those of other SUSY breaking 
parameters. U (1) b-l extension of the MSSM might provide the origin of large Bn- It would 
be interesting if such a large value oi Bn can be naturally possible. 
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Appendix A: CP violation and phase convention 

Here, we discuss the CP violation of SUSY type I seesaw model and identify the inde- 
pendent phases by choosing a phase convention. One can assign the R charge to the Higgs 
superfields Hi and H2, and 1 to the lepton superfields L and N"^. Under the R transfor- 
mation and the phase redefinition of the superfields L, Hi, H2, N'^, the super potential is 
transformed as, 



W ^ e^*^« y-Y, exp(z(^;v= + 0^ + e2))N'L ■ H2 - ^ exp(2«0jvO^'^'+ 

^iexp{t{ei + 02- 2eR))Hi ■ H2)) . (Al) 



Therefore one can remove the phases of the parameters YJ^, /i, Mr in W by choosing the 
phases of the superfields as follows, 

On'^ = argMfl, 

61 + 62- 26r = -arg/i, 

6L + 62 + 6N^ = -aTgY,, (A2) 

The trilinear couplings of the soft breaking terms transform in the same way as the super 
potential, so one can not remove those phases. For the soft breaking parameters of the 
bilinear form, one can take one of them to be real. We then rotate the phase of away 
by choosing the phase parameter of R transformation as follows, 

6R = ~aTg{BN). (A3) 

In Eg. (lA2p . we still have the freedom of choosing the phase of 62- Here, we choose the phase 
62 so that vacuum expectation value of H2 becomes real 



^2 = -arg(t;2). (A4) 



To summarize we choose the phases as. 



6'Arc = argMi? 
61 = — arg fi + arg(t'2) — arg B 



N 



6l = aTg{v2) + ^ arg Mr - arg Y^. (A5) 
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With this phase convention, the soft breaking terms are written as, 



■'soft 



AA\YAN*e 



+ h.c. 



+ 2\n\\B\Re.{e'^^^^^Hi ■ H2) - 



R\ 



Bn\N*N* 



- ml\L\'-m%\N\\ 
and two independent irremovable CP violating phases are presented as, 

B = \B\e''''^^, 



(A6) 



(A7) 



Appendix B: derivation of the effective potential 

In this appendix, wc derive the effective potential of Higgs fields in SUSY type 1 seesaw 
model. The contribution to the effective potential for Higgs fields from the loops mediated 
by neutrino superfields is written as, 

Q^-'^d'^k 1 



(lndet(Ms2 - k^) - lndet(MF , 



(Bl) 



{2nyi 2 

where Mp is the mass matrix of one of the neutrino sector and is the 4 by 4 mass-squared 
matrix of sneutrino sector given by 

f (mj+ml) AlniD IMruid] \ 





A^niD 



\mnMR\ \MR\^ + ml {B^MrI 
AlrriD \B*t,MR\ \MR\'' + m\) 



(B2) 



where 







rriD 




A. 








A. 





(B3) 



The effects of CP violation appear through the parameter A^,. We compute the following 
quantity. 



lndet(Ms^ - k^) = TrLn(Ms^ - k^). 
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(B4) 



To compute the scalar contribution, we diagonalize M| approximately and treat A term as 
perturbation. We first split as 



where 



f (m| + m%) 







(ml + ml) IMruid] 

IuIdMrI \ivir 



MR\^ + m% + ml, 



N 







B%Mr\ 



IMrTUdI 



\BnMr\ 
\MR\^ + ml + mlJ 



and 



/ 


A,,mn 
V AlrriD 



AlruD \ 
A^niD 


J 



(B5) 



(B6) 



(B7) 



One can find the orthogonal matrix O which diagonahzes Mg. Using this matrix, MJ^ is 
transformed as 



OMs'O' = diag(ml, m|, m^, mi) + OA^O' . 



(B8) 



Here, mi,m2 are the mass of lighter sneutrinos and m^^mi are those of heavier sneutrinos 
given by 



Ml 


+ m\ 


+ 2m|) 


+ BnMr 


+ m| 


Ml 


+ m\ 


2 


-BnMr 


+ 

Li 


_ Ml 


+ ml 


2 

+ 2m|) 


-BnMr 


+ mp 

Li 


_ Ml 




2 

+ 2m|5 


+ BnMr 


+ m| 






2 







{Ml + m\- B^Mr - m|)2 + AmlMl, 



+ \\J{Ml + ml + BnMr - m\y + AmlMl. 

(B9) 



These sneutrino masses should be compared with the neutrino masses written as. 



H 



rrir = 



R 



Ml 



2 VM^ + 4m^M| 



(BIO) 
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Using Eq.( ]Bl|) and the mass eigenvalues, one can find effective potential as follows, 



where 



= ; / ^^7^^ I > ;ios(™? - *') - 21og(m-J, - e) - 21og(mi - e) 




^, i=l 



(B12) 



where C[/y = 7 — 7 + log 47r and Q is the renormalization scale. The renormalization point 
dependent finite part of the effective potential V^^ is given as, 

"'"'(O^) - ^ [t '"?('°^ f - 5) - 2<(log f - 2) - 2»i(log f - f )) (B13) 

We note that V^^^ depends on the Higgs vacuum expectation value through mj^ where = 
To obtain the contribution to the Higgs mass term m\j^H\H2-, one can differentiate 
the effective potential with respect to m|,, while keeping the terms which remain non zero 
in large limit of M/j, 



2 



1 M| ^ ^dml ^dml ^dmrr 

1 / \ 1 

^ T^((l°gQf)K + + 2^^j-Y^(C^v + l)K + m^), (B14) 

The terms which are proportional to the derivative of the lighter mass also vanish in large 
limit of Mn, because ml ~ ml — wi| ~ and the derivatives with respect to mj^ are 



suppressed as ^ and jp-, respectively. From Eq. ( ]B14p . one can read off the coefficient of 
the Higgs mass term H\H2. The contribution to the Higgs mass term including the counter 



term is given as 



r2 



^,{HlH,) [\og -^{ml + ml) + 2ml ) (B15) 
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where the counter term is given as, 

Y 2 



{Cuv + l){m)+ml)HlH2. 



(B16) 



Next we compute the corrections to V^^*^-* due to the terms up to the second order 
of A^, because they give the non- vanishing contribution to the effective potential in large 
limit of Mr. To compute the corrections, one needs to derive the orthogonal matrix O in 
Eq.( ]B6l) . To diagonalize Mq follows two steps. First, we diagonalize Mq with the help of 
orthogonal matrices Ol and Oh as follows. 




01 







OS 







m| + mi) 





TUdMr 





tudMr Ml + m% + ml- B^Mr 

\ uidMr 



mDMR \ 



M^ + m% + ml + BnMr ) 



(BIT) 



where Ol and On are given as 



1 1 
-1 1 



(B18) 



We note the degenerate diagonal masses of the heavy sneutrinos are split after the rotation. 
The mass squared matrix M'^ has the separated two by two parts as sub-matrices. Each of 
them has the form of the seesaw type. Thus, the mass matrix Mq can be diagonalized as. 



1 








\ 









m| 



\ ^ m\j 

( cose+ 

COS 6^ — sin 6^ 
sin 6_ cos 6_ 
\ sin ^+0 



- sin 6*+ \ 



cos 9-1- 



( cose+ sine+\ 

cos^„ sin^^_ 

— sin 6_ cos 6_ 
\ - sin 61+ cos 6+ ) 



.(B19) 
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Then the orthogonal matrix O is given as 

/ cos^+ 
cos^_ — sin^_ 



X 




sm( 


sine_ cose_ 
V sin 6'+ cos 9+ J 

Using the above form of orthogonal matrix O, OA^O^ is given as 

/-sin 2^+ cos2^+\ 

sin2^_ -cos 2^^ 

-cos2^_ -sin2^_ 

\ cos 29+ sin 29+ J 

( sin(^_ + e+) -cos(e_ + ^+) 

-sin^_+^+ cos(^_ + 

cos(^_ + ^+) sin(^_ + 

V -cos(^_ + ^+) -sin(^_ + ^+) 



(B20) 



OAaO^ = mDRe(i^) 



+ 



\ 



) 

(B21) 



We then obtain the corrections to the effective potential at the first order of given as 



2 l^^i {2nyi mf-k^ J 



RcA^niD 



327r2 



X 



(2 2 
ml sin 29+{Cuv + 1 - In ^) - sin 29+ {Cuv + 1 - In ^) 



m 



- ml sin 29- {Cuv + 1 - In — |) + m2 sin 29^ {Cuv + 1 - In 



mr, 



(B22) 



Q'' ' g2. 

Now, let us show how the divergences are canceled so that the correction is finite. To do 
this, we use the relation 



(777.4 ~ "^1) sin 29+ — {m^ — 7712) sin 29- 



(B23) 
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Then, the corrections to the effective potential become 



SV^ = ^^^^^ ml sin 26^ In ^ - sin 26. In 



777-3 20- + m4 sin 2^+ m| 

^ 2 



327r2 

~ ^ (i?e.(A,5^)i/ii/2 - ^iB^Re.iH, ■ H^)) , (B24) 

where we have used the relation which is valid in large hmit of Mr, ^± ~ ^ and m| — m| = 
2MjiB]sr. The correction at the second order of A^^ term is given as 



The term which is not suppressed by is given as. 



(B26) 



The divergences are canceled by adding the counter term, 

K^') = ^{Cuv + 1) [\A,\^hI ■ H2 - 2Re{A,fiH, ■ H^) + //'//j • . (B27) 
The effective potential at one loop level is finally written as 

a^^(g2)) h\H, + + ml^{Q^)) hIh^ - 2Re.{B{Q^)fj.H, ■ H^) 

gf K + < + 1^' I) + 2< + 2Re.(A,5^)J i/ti/2 







+ 








+ 




167r2 




2Re. 



where >Cd is the D-term contribution. To complete the renormalization of the effective 
potential, we consider the relation between the renormalized mass parameters and the bare 
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ones. We first note that the bihnear part of the Higgs sector including the tree and the 
counter terms in the present model can be derived from the following Lagrangian, 

£ = ZiHIHiId + Z2HIH2\d + f^Hi ■ H2\f + h.c. 

+ 2Re. ((5(g2) + SB)fiHi ■ H2) . (B28) 
After integrating out F terms of the superfields, one obtains, 
£ = Zid^Hld^Hi + Z2d^Hld^H2 



H\Hi - ^-^hIH2 + 2Re. ((^(Q^) + 5B)iiHi ■ H2) 



Z2 Z\ 

- + 5ml_MH, - imjj^iQ') + 5mljHlH2. (B29) 

We define bare superfields and bare parameter fi as Hf = y/^iHi, {i = 1, 2) and 
lii^\fZ{\fZ2 = /i, respectively. One can write the Lagrangian in terms of the bare fields 

as, 



^uT^u _ my-^ ^ ■ ^ ""^ H^^Hl (B30) 
Zi Z2 



Then one can define the bare mass parameters as 



Bo = B{Q^) + 5B. (B31) 



Eq. (IB29| ) leads to the following counter terms for the bilinear parts of the Higgs potential, 

K = {5ml^ + {Z^^ - l)\fi\'')HlH, + + (^1"' " IM')hIH2 

- 2Re{5B^jiHi- H2). (B32) 

Comparing Vc with the sum of the counter terms V^^^"* + vi^"* given by 

+ = ^,{Cuv + MAA' + ml + m^^HlH2 

+ ^2(Cuv + l)^^'HlH, 

-2 -^{Cuv + l)Re{A,^^H^ ■ H2), (B33) 
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we obtain the following relations, 

Sm%^ + (Z^' - l))"' = ^2iCuv + l)(|^.r + + ml), 

5B = -Ji^(Cuv + l)A.. (B34) 
Using the results of the wave function renormalization, 

Zi = 1, 

iS^^^^' ^^^^^ 



we obtain 



Sml^ = ^/^^ (B36) 



^""^^ = iS^^''^'^'' + + + !)• (B37) 



Finally, we find the following relations between the renormalized parameters and the bare 

ones, 

2 //o2\ 2 ^ (/ 2 

m^^(Q ) = moH, - Y^l^ , 

i?(g2) = 5o-^A,(C[/y + l), 

l2 



i^{Q') = /^oV^2. (B38) 
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